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ABSTRACT 

We study galaxy evolution from z = ltoz = 0asa function of velocity dispersion 
ct for galaxies with a > 95 km s _1 based on the measured and Monte Carlo realised 
local velocity dispersion functions (VDFs) of galaxies and the revised statistical prop- 
erties of 30 strongly-lensed sources from the Cosmic Lens All-Sky Survey (CLASS), the 
PMN-NVSS Extragalactic Lens Survey (PANELS) and the Hubble Space Telescope 
Snapshot survey. We assume that the total (luminous plus dark) mass profile of a 
galaxy is isothermal in the optical region for ^ z ^ 1 as suggested by mass modelling 
of lensing galaxies. This study is the first to investigate the evolution of the VDF shape 
as well as the overall number density. It is also the first to study the evolution of the 
total and the late- type VDFs in addition to the early- type VDF. For the evolutionary 
behaviours of the VDFs we find that: (1) the number density of massive (mostly early- 
type) galaxies with a > 200 km s _1 evolves differentially in the way that the number 
density evolution is greater at a higher velocity dispersion; (2) the number density 
of intermediate and low mass early-type galaxies (95 km s _1 < a < 200 km s _1 ) is 
nearly constant; (3) the late-type VDF transformed from the Monte Carlo realised 
circular velocity function is consistent with no evolution in its shape or integrated 
number density consistent with galaxy survey results. These evolutionary behaviours 
of the VDFs are strikingly similar to those of the dark halo mass function (DMF) 
from N-body simulations and the stellar mass function (SMF) predicted by recent 
semi-analytic models of galaxy formation under the current ACDM hierarchical struc- 
ture formation paradigm. Interestingly, the VDF evolutions appear to be qualitatively 
different from "stellar mass-downsizing" evolutions obtained by many galaxy surveys. 
The coevolution of the DMF, the VDF and the SMF is investigated in quantitative 
detail based on up-to-date theoretical and observational results in a following paper. 
We consider several possible systematic errors for the lensing analysis and find that 
they are not likely to alter the conclusions. 
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1 INTRODUCTION 

The observationally-derived statistical properties of galax- 
ies provide key constraints on models of galaxy formation 
and evolution. The statistical properties of galaxies include 
the luminosity functions (LFs), the stellar mass functions 
(SMFs) and the velocity functions (VFs). In the current 
Lambda cold dark matter (ACDM) hierarchical structure 
formation picture CDM haloes are formed due to gravita- 
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tional instabilities and evolve hierarchically th rough merging 
fe.g. lWhite fc Rees!ll978l : lLacev fc Coleill993l ). Baryons set- 
tle in the dark matter halo potential wells and undergo dis- 
sipational radiative processes including star formations and 
evolutions, supernovae explosions, AGN activities and their 
feedback, resulting in visible galaxies (in the central parts 
of the haloes). Visible galaxies may further merge (and be- 
come morphologically transformed) resulting in evolutions 
of galaxy populations. The various statistical functions of 
galaxies are the end results of complex processes involving 
gravitational physics, baryonic physics and merging. 
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The LFs of galaxies have traditio nally been mea- 
sured most extensivel y and reliably (e.g. Cole et all 200 ll; 
Kochanek et al.l l200ll: iNqrberg et all 120021 : iBlanton et all 
20031 ; ICroton et al.l I2005T ) and consequently have pro- 
vided vital constraints on models of galaxy forma- 
tion and evolution. More recently, the SMFs have been 



the focus of many observational studi e s (e.g. Cole et all 



200 ll; iBell et al.l 120031; iDrorv et all 12004 iBorch et al.1 12006 ; 
Bundv et al.ll200d; ICimatti et al.ll2006l; iFontana et al lbood; 
Pozzetti et al. |2007|~ Conselice et all 120071 ; IScarlata et all 
20071 ; iMarchesini et alj |2009l because they can provide more 



robust tests of the underlying ACDM theory. Recent semi- 
analytic models o f galaxy formation and evolution (e.g. 
Croton et af] 120061: I Bower et alllioOr] ; iKitzbichler fc White! 



20071 ; IStringer et all I2009F appear to match the observed 



local LFs/SMFs owing in large part to the suitable incorpo- 
ration of baryonic physics such as AGN feedback. 

However, the evolutions of the LFs/SMFs from a high 
redshift Universe are not understood as well. On the ob- 
servational side, deep surveys of galaxies have been used 
to derive the evolutions for the type-specific populations 
as well as the total population, particular attention paid 
to photometrically red galaxies or morphologically early- 
type galaxies which occupy the most massive part of the 
galaxy population. Various existing observational results on 



survev results (e.g. ICimatti et al.l 20061: Bundv et al.1 2006: 


Fontana et al.1 l2006t iPozzetti et al.1 12007 


; IConselice et al. 


2007 




Scar lata et al.l 12007; Brown et al.l 


2007; Cool et al. 


2008 




Marchcsini et al.l 20091) suggest little evolution of 



massive (super-L*) early-type galaxies often accompanying 
"mass-downsizing" behaviour in which more massive galax- 
ies are to be in place at an earlier time so that the num- 
ber density of most massive early-type galaxies evolves least 
while the number density of typical ~ L* galaxies may 
evolve significantly over cosmic time. On the other hand, 
there are also galaxy survey results that do not partic- 



many support, mass-aownsizmg evolutions ^e.g. r>en er au 
2004; lilbert et al.ll2006l; iFaber et al.ll2007l; lLotz et al.1 2008; 


Ilbert et al.1 120091). On 
model predictions on 
(e.g. IBower et al.1 12006 


the theoretical side, semi-analytic 
the LF/SMF evolutions scatter 
; iMenci et al.1 120061; iMonaco et al.l 


20061; de Lucia et all 


2006; Kitzbichler & White 20071; 


Almeida et al. 


2008:; ICattaneo et al.1 20081; Stringer et al.1 



2003 ; iFontanot et al.l 12003 ) and are not in good agree- 
ment with existing observational constraints on the evolu- 
tions of the SMFs, particularly for the type-specific galaxy 
populations. Semi-analytic models do not predict mass- 
downsizing evolutions of the SMF and according to more 
recent studies the mismatch between model predictions and 
observed SMFs is more severe for less massive galaxies with 
M star s < 10 n M w (e.g. ICattaneo et al.ll2008l ; IStringer et all 



2009!; IFontanot et al.ll2009l ) 

Given the varying results on galaxy evolutions both 
observationally and theoretically, it would be of great 
value to constrain galaxy evolutions through other inde- 
pendent methods based on different data sets. The sta- 
tistical analysis of strong lensing galaxies in conjunction 
with well-determined local velocity dispersion functions 
(VDFs) of galaxies provides such an independent method. 
The VDFs represent the statistical properties of the dy- 
namics of galaxies that have been modified by baryonic 



physics from the CDM haloes. The velocities of particles 
and stars are boosted and the density profiles become 
steeper. These dynamical features can be seen by compar- 
ing the dynamical properties of A-body simulated haloes 
and the observationally-derived kinematical and dynami- 
cal properties of visible galaxies. N-body simulations show 
that pure CDM haloes follow NFW ( or NFW-like ) den- 
sity profiles (e.g. INavarro et al. _ 1997 : iMoore et al.l 1 19991 ; 
IJing fe Sutcl 120021 ; INavarro et all I200I ) and hence rising 
velocity profiles in the inner parts whereas the kinemat- 
ical and dynamical properties of galaxi es, derived from 
such methods as stel lar dynamics (e.g. iRix et aH 1 19971: 
Cappellari et al.l 120061). rotation curves (e. g. Rubin et al.l 
19851 : iPersic et al.l 1 19961; Salucci et al.l 2007 | ). gravitational 



lensing (e.g. Rusin fc Kochanek! 120051 ; IrCoopmans et al.1 
l2006l ; lG"avazzi et al.ll2007h . show that baryon settled galaxies 
follow isothermal (or isothermal-like) density profiles imply- 
ing that velocities are boosted and mass distributions be- 
come more concentrated in the inner regions. An important 
observable quantity that characterises the dynamical state 
of the galaxy is the line-of-sight velocity dispersion a es- 
pecially for rotation-free early-type galaxies. For late-type 
galaxies, rotations are important in the disk. Nevertheless, 
the bulge and the surrounding halo of the late-type galaxy 
may be regarded as an ellipsoidal system that is similar to 
early-type galaxies. Furthermore, given that most late-type 
systems are isothermal(-like) as inferred from nearly flat ro- 
tation curves, we may relate the circular rotation speed in 
the disk v c to the velocity dispersion in the inner halo via 
v c = V2a. 

All three statistical functions (i.e. LF, SMF, and VDF) 
are the consequences of baryonic physics from pure CDM 
haloes. However, unlike the other functions the VDF has 
only to do with the dynamical effects of baryonic physics 
separated from many other complex effects involving radia- 
tive processes. Furthermore, given that galaxy mass profiles 
are likely to be isothermal(-like) up to significant fractions 
of the dark halo virial radii as stressed above, we expect a 
good correlation between the central velocity dispersion and 
the halo virial mass. Hence the velocity dispersion function 
(VDF) will be another independent and powerful constraint 
on galaxy formation and evolution. One advantage of the 
VDF as a cosmological probe is that it can test the un- 
derlying ACDM structure formation theory more robustly 
because it is linked more intimately to the theoretical halo 
mass function. 

In this paper we use a range of local VDFs. The local 
VDFs include n ot only the SPSS VDF of early-type galax- 
ies measured by I Choi et al.l l|2007h but also the type-specific 
and total VDFs obtained through a Monte Carlo method 
based on the galaxy LFs from the SDSS and intrinsic corre- 
lations between luminosity and velocity dispersion or circu- 
lar rotation speed. We then constrain the evolutions of the 
VDFs through the statistical analysis of galactic-scale strong 
gravitational lensing. Strong lensing statistics is a powerful 
means to constrain the evolution of the VDF because strong 
lensing probability is proportional to comoving number den- 
sity times a 4 (the image separation is proportional to a 2 ) 
and lensing galaxies are distributed over a range of redshifts 
(up to z > 1) and virtually over all area of sky. The latter 
property is particularly important since it means that unlike 
deep galaxy surveys strong lensing data are free from cosmic 
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variance. For strong lensing statistics we use a total of 30 
lenses fro m the radio-selected Cosmic Lens All -Sky Survey 
(CLASS; iMvers et al.ll2003l ; iBrowne et al]l2003l) and PMN- 
NVS S Extragalactic Lens Survey (PANELS; IWinn et all 
l200lT ) and the o ptically-selected H ubble Space Telescope 
Snapshot survey l)Maoz et al.l Il993h . as these surveys have 
well-defined selection functions and cover about the same 
image separation range (i.e., > 0".3) and similar redshift 
ranges for source quasars El However, only the CLASS sta - 
tistical sample of 13 lenses (|f3rowne et al.ll2003l ; IChadl2003l ) 
and the Snapshot sample of 4 lenses are used for the statis- 
tics of absolute lensing while the rest of lenses are used only 
for the distributions of image separations (and lens redshifts 
if available). 

This paper is organised as follows. In §2 we derive 
through a Monte Carlo method the local VDFs for both the 
type-specific galaxy populations and the total population. 
In §3 we describe the strong lensing data and the theoreti- 
cal model for statistical analysis of strong lensing. In §4 we 
constrain (non-evolving) intermediate redshift VDFs using 
strong lensing data and then compare the intermediate red- 
shift VDFs with local VDFs. This simple analysis captures 
the trends of the evolutions. Next we study in detail the 
evolution of the VDF through a parametric approach (§5). 
The results that come out are striking. The number den- 
sity evolution as a function of a is differential: the number 
density of galaxies of typical velocity dispersion and lower 
evolves little since z — 1 but that of larger velocity disper- 
sion evolves significantly showing "velocity upsizing (hier- 
archical)" behaviour of galaxy evolution. In the discussion 
section we first consider possible systematic errors (§6.1) 
and then compare with previous results of strong lensing on 
galaxy evolutions (§6.2) and the evolutionary behaviours of 
the DMF and the SMF from the literature (§6.3). A com- 
prehensive analysis of the coevolution of the DMF, the VDF 
and the SMF is carried out in a companion paper l|Chael 
2010). We give the conclusions in §7. Throughout we as- 
sume a ACDM cosmology with (S7 m o,r2Ao) = (0.25,0.75) 
and H = l OO/ikms" 1 Mpc" 1 co nsistent with the WMAP 
5 year data ijDunklev et al.l [20091 ). When parameter h does 
not appear explicitly, h — 0.7 is assumed. 



2 VELOCITY DISPERSION FUNCTIONS OF 
GALAXIES AT Z w 

2.1 Statistical functions of galaxies 

The Schechter luminosity function (LF) 0l, the differen- 
tial comoving number density as a function of luminosity, is 
given by 



dn = (f>-L{L)dL = 



L" 



exp 



L* 



dL 



(1) 



It has been measured extensively at various wavebands and 
rece nt observations have produced reliable results at optical 



(e.g. iNorberg et al.l :2002: Bla nton et al . 2003; C roton et al 



l2005h and infrared (e.g. ICole et all l200ll ; iKochanek eTal 



1 SDSS quasar lens search l llnada et al ] |2008h is limited to image 
separation range of 1" and source redshift range of 0.6 < z < 
2.2. 



l200ll) wavebands owing particularly to large surveys such as 
the Sloan Digital Sky Survey (SDSS), the Two Degree Field 
Galaxy Redshift Survey (2dF) and the Two Micron All Sky 
Survey (2MASS). 

Given the observed power-law correlations between lu- 
minosity and internal velocities of galaxies such as the Faber- 
Jackson llFaber fc Jacksonl Il976h relation and the Tully- 
Fisher (|Tullv fc Fisher! Il977l ~ relation, one may expect a 
modified for m of the Sche c hter f unction for a velocity func- 
tion. Indeed, ISheth et al] (|2003h find that the distribution 
of the velocity dispersions of early-type galaxies from the 
SDSS is well fitted by the VDF of the form 



dn = <j>{a)d<7 = <j>* I — J exp 



o_\P 
a" 



f3 da 



(2) 



or equivalently 

<!>(V)dV = ri0 aiV - yt) exp [-lO^-^'l i^dV, (3) 

L J 1 (a/p) 

where V is the logarithmic variable given by 



V = lo Sio G - 



(4) 

Notice that equation @ along with equation {TJ imply the 
following correlations: 



Ql = a/0-1 and 4>t = cj>* /T(u/p), 



(•>) 



where <j>* is the integrated comoving number density of 
galaxies, and an 'effective' power-law correlation between 
luminosity (L) and velocity dispersion (er) of 

L_ 

17 = 



(6) 



However, directly measured relations of L = L(a) and 
a — a(L) give different values for the power-law slope j3 
compared with that in equation (|6} because the scatter in 
the correlation ac ts in different ways for the different fits 
jSheth et all 120031 ). Furthermore, the value of (3 varies as 
a function of luminosity as can be seen in a magnitude- 
ve locity dispersion plane of galaxies [see Fig. [1] or Fig. 4b 
of IChoi et all (|2007h ]. This is why it is not reliable to turn 
a LF of galaxies to a VDF simply using equation (J5J) and 
a (mean) Faber-Jackson or Tully-Fisher relation although 
this was often practised in the past. One must fully take 
into account the intrinsic scatter in the magnitude-velocity 
dispersion plane. This is done below. 



2.2 VDF of early- type galaxies 

IChoi et ail (|2007l ) have measured the VDF for local early- 
type galaxies based on a large sample of galaxies from 
the SDSS DR5 data set, employing a galaxy classifi- 
cation scheme that closely matches visual classifications 
llPark fc Choil 120051 ). The early-type VDF measured by 
IChoi et al] (|2007l ) will be referred to as 'E0' VDF and can 
be found in Table [1] Notice, however, that galaxy number 
counts start to become incomplete at lo w velocity disper - 
sions for a given galaxy sample used by IChoi et aD (|2007l ) 
because the sample is n ot limited b y veloc ity dispersion but 
by absolute magnitude. IChoi et all (|2007l ) discuss this kind 
of problem for a < 150 km s _1 and obtain galaxy counts 
down to a ~ 70 km s _1 using a series of samples of differ- 
ent volume and magnitude limits. The measurement of the 
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Figure 1. Monte Carlo realised magnitude- velocity dispersion (M r — a) relations. Only 5% of the gener ated galaxie s are s hown. The 
magnitude is the SDSS r-band absolute magnitude. The relation for the early-type population is from the Choi et al measured 
relatio n while that for the late-type population is inferred from the relation between M r and the rotational circular speed bv lPizaeno et all 
l|2007h assuming a singular isothermal sphere for the total mass distribution of a galaxy. The relation for all galaxies is the combination 
of the early- type and the late- type relations. 



LF is free of this problem at low luminosities down to the 
absolute magnitude limit of a given sample. Hence galaxy 
counts are more reliable at low luminosities than low velocity 
dispersions. 

It is therefore desirable to use galaxy counts at low lu- 
minosities in deriving the VDF for the full range of velocity 
dispersions. We use a Monte Carlo method as follows. We 
generate galaxies within a comoving volume brighter than 
an absolute magnitude limit according to a measured lumi- 
nosity function. The next step is to assign a velocity disper- 
sion to each galaxy of a given absolute magnitude. To do 
this we use the 2-dimensional distribution of galaxies in the 
absolute magnitude-velocity dispersion (M — a) plane. We 
assume that the most likely value of V = log 10 a varies as 
a function of M r (absolute mag nitude in SDSS r band) ac- 
cording to the lChoi et al.l (|2007l ) Fig. 4b. We assume that V 
follows a normal distribution for a given M r with a varying 
dispersion of 



s{M r ) = a(M r - M*) + s* , 



(7) 



as suggested by the lChoi et al.l l|2007h Fig. 4b. Once all mock 
galaxies are assigned velocity dispersions, we finally derive 
the VDF by counting mock galaxies as a functio n of V. 

W e cons ider the SPSS ea rly-type LF by IChoi et ail 
(|2007l ). The IChoi et al.l (|2007h LF is the opt i mal choice 
for the following two reasons. First, IChoi et"al] (|2007l ) clas- 
sify galaxies using a class ification method t hat matches 
closely visual clas sification (|Park fc Choi|[2005l ) . Second, the 
IChoi et all (120071 ) LF is based on the same SDSS magnitude 
system that is used for the magnitude-velocity dispersion 
relation, so that any error arising from photometric trans- 
formations (e.g., from the 2dF magnitude) is avoided. For 
this LF we derive the early-type VDF using the Monte Carlo 
method described above adopting M* — 51og 10 h = —20.23, 
s* = 0.085 and a = 0.025 where s* and a are chosen so 
that the simulated M r — V relation mimics that observed 



by IChoi et"al1 (|2007r i with a mean dispersion of ~ 0.11 
for —22 < M r < —17.5. The simulation comoving vol- 
ume is 10 7 h~ z Mpc 3 and the faint magnitude limit is 
M r — 51og 10 h — — 16. This simulation generates more than 
100,000 ear ly-type galaxies. T he simulated M r — V rela- 
tion for the IChoi et all (|2007f l early-type LF can be found 
in Fig. [1] Notice that th is simulated plane is very similar 
to the IChoi et al.l (|2007r i Fig. 4b. The derived early-type 
VDF will be referred to as 'El' VDF and can be found 
in Fig. [2] and table [1] Notice that the VDF derived from 
the IChoi et al.l (|2007r) early-type LF closely matches the 
IChoi et all (|2007V ) measured VDF for a > 150 km s" 1 but 



shows significant difference at lower a. This could mirror 
the inco mpleteness in gala xy counts at low velocity disper- 
sions by IChoi et al.l (|2007T> . Becau se of this possibility we 
consider both the Choi et all l|2007h measured VDF and the 
Monte Carlo realised VDF. 



2.3 VDF of late- type galaxies 

There does not yet exist any kind of directly measured ve- 
locity functions for the late-type population of galaxies. For 
the late-type population, the internal motions in the optical 
regions of galaxies are observed as rotations since the (ran- 
dom) motions of dark matter particles are not observed. Un- 
fortunately, it is not an easy observational task to measure 
the rotation of a galaxy and consequently rotations of galax- 
ies have not been measured for all galaxies in a volume or 
magnitude limited galaxy sample. However, one may obtain 
a velocity function for the late-type population from the LF 
through the Tully-Fisher relation between rotational speed 
and luminosity that is usually a tighter correlation than the 
Faber- Jackson relation o f the e arly-type population. For ex - 
ample [Gonzalez et all il2000h iKochanek fc Whltel i|200ll ). 
IChad i|2003h and lSheth et all l|2003h inferred circular veloc- 
ity functions (CVFs) of late- type galaxies combining various 
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Table 1. Local statistical functions of galaxies based on SDSS data 



Type 


ID 


a 







<f>* 

(h 3 Mpc" 3 ) 


(h 3 Mpc" 3 ) 


M * — 5 log 10 h ct* 

( km s — 1 ) 


Reference 




Early-type 


VDF 


EO 


2.32 


2.67 


-0.13 


0.008 


0.0073 


161. 


Choi ct al. 


(2007) 


LF 








-0.527 




0.0071 


-20.23 


Choi ct al. 


(2007') 


VDF 


El 


0.85 


3.72 


-0.77 


0.013 


0.0033 


217. 


This work 




Late-type 


LF 








-0.897 




0.012 


-20.22" 


Choi ct al. 


(2007) 


VDF 


LI 


0.69 


2.10 


-0.67 


0.066 


0.024 


91.5 


This work 












a' 




All 

e 








VDF 
VDF 


AO 
Al 


0.69 
0.69 


2.01 
2.05 


6.61 
7.69 


0.074 
0.078 


0.044 
0.035 


100. 
102. 


This work 
This work 





Notes. a An internal extinction of A r ~ 0.1 [extrapolated from Fig. 11 of IPizagno et aT I ll2007l) l is assumed and corrected for the 

IChoi et ail (|2007t l b/a > 0.6 late-type galaxies. 



late-type LFs and Tully-Fisher relations. It is also possible 
to convert the CVF to the VDF via a = n c /\/2 assuming 
the singular isothermal sphere model for the total mass dis- 
tribution of a galaxy0 In this work we essentially follow this 
approach to obtain estimates of the VDF of late-type galax- 
ies. However, we fully take into account the intrinsic scatter 
in the Tully-Fisher relation using the Monte Carlo method 
described in §2.2. 

The Tully-Fisher relation has been published by numer- 
ous authors. The relation is usually expressed as L/L* = 
(vc/vc) 13 or its inverse relation v c /v^ = (L/L*) 1 ^. The slope 
P ranges from ~ 2.5 - 3.5 in most of the published res ults 
in optical and infr ared bands (Se e IPizagno et all 120071 and 
references therein) . IPizagno et al.l (|2007 ) derive the relation 
L = L(v c ) and its inverse relation v c = v c (L) for 162 SDSS 
ga laxies in the SPSS g. r, i and z bands. The galaxies used 
bv lPizagno et all (|2007l ) have apparent minor-to- major axis 
ratio b/a ^ 0.6 so that internal extinctions are measured and 
corrected for reliably. The measurements of the TF relation 
and its inverse relation give different slopes because of the 
scatter (as for the FJ relation). Notice that the inverse rela- 
tion v c /v* — (L/L*) 1 / 13 is more applicable in deriving a ve- 
locity function from the LF because the velo city function is 
estima ted from the LF. For the SDSS r-band IPizagno et al.l 
<|2007| ) fmdlog 10 (v c ) = (-0.135±0.006)(M r -M r *) + (2.210± 
0.006) with M* — 51og 10 h — —20.332 and a dispersion of 
0.063 ± 0.005. The relation between V(= log 10 a) and M r 
then follows from log 10 (w c ) = V + \og 1Q y/2 for the SIS model 
of the galaxy. 

Using the Mon te Carlo method based on the above 
IPizagno et~ai1 l|2007t ) M r — V relation th e VDF of late-typ e 
galaxies is derived for the late- type LF by I Choi et al.l (2007). 
As in §2.2 the simulation comoving volume is 10 7 h~ 3 Mpc 3 
and the faint magnitude limit is M r — 51og 10 h = — 16. By 



2 Dark matter haloes are, in general, triaxial in shape and have 
density profiles that deviate from the isothermal profile. Never- 
theless, the deviations of shapes and profiles from the SIS may 
have cancelling effects statistically so that the averaged property 
of a galaxy might not deviate too much from the SIS. 



this simulation more than 300,000 late- type galaxies are gen- 
erated. The derived VDF will be referred to as 'LI' VDF 
and can be found in Fig. [2] and Table [T] Fig. [2] a lso shows a 
VDF derived assuming no intrinsic scatter for thelCh oi et al.l 
l|2007l) late-type LF to demonstrate the effect of the intrinsic 
scatter in the Tully-Fisher relati on. The simulated M r — V 
relation for the lChoi et al.l l|2007h late-type LF can be found 
in Fig. [1] 



2.4 VDF of all galaxies 

In §2.2 and §2.3 the VDFs are derived for morphologically- 
typed galaxies. One may well expect that a well-defined 
VDF exist for the entire population of galaxies. Here we con- 
sider this possibility. The total VDF is motivated in several 
respects. First of all, the total VDF can be compared with 
theoretical predictions more straightforwardly than type- 
specific VDFs just as is the total LF because galaxy clas- 
sifications are not required. Secondly, the evolution of the 
total VDF is of interest in its own right as it represents 
the collective evolution of all galaxies including morphologi- 
cal transformations through merging. Thirdly, strong lensing 
computation is simplified if the total VDF is used because 
one does not require knowledge of the lensing galaxy type 
which is not well determined by observations for some sys- 
tems. Lastly, the observed evolutions of the total LF/SMF 
may be used for or compared with the evolution of the total 
VDF. 

It is straightforward to derive the total VDF using the 
Monte Carlo method used above. Early-type and late-type 
galaxies are separately generated and then combined in the 
M r — V plane. For the comoving volume of 10 7 h~ 3 Mpc 3 
and the faint magnitude limit of M r — 5 log 10 h = — 16 over 
400,000 velocity dispersions of galaxies are use d to derive 
the to tal VDF. Again we consider the LFs by IChoi et"ai1 
(2007). We also con s ider t he early- type VDF directly mea- 
sured bv lChoi et ail |2007t) . The derived total VDFs can be 
found in Fig. [5] The 'AO' VDF corresponds to the combina- 
tion of the E0 VDF and the LI VDF while the 'Al' VDF 
corresponds to that of the El and the LI. Notice that the 
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7=LOG 10 (cr/km s" 1 ) 

Figure 2. (Upper panel) The local velocity dispersion functions 
of galaxies based on SDSS data (see also Table [T}. The red, blue 
and black curves represent respectively the VDFs for the early- 
type, the late- type and all populations. The dashed curves are 
Monte Carlo realised VDFs. The r e d soli d curve is the early-type 
EO VDF measured bv lChoi et"al] ll2007t) . The black solid curve 
is the total VDF resulting from the combination of the EO VDF 
and the Monte Carlo realised late-type VDF. The blue dotted 
curve is the late-type VDF that would be obtained ignoring the 
intrinsic scatter in the magnitude- velocity relation. (Lower panel) 
The behaviours of 4>(V) X cr 4 where c/>(V) is the differential number 
density of galaxies and tr 4 is proportional to the lensing cross 
section of each galaxy (assuming the isothermal galaxy profile). 
These curves show the relative lensing efficiencies of the early- 
type and the late-type galaxy populations as a function of velocity 
dispersion. 



simulated data points for the total population of galaxies 
are not fitted well by the VDF of equation ([2]) in contrast to 
the morphologically-typed galaxy populations. This might 
be the true nature of the total VDF. However, it could be 
due to the errors in the adopted correlations between lumi- 
nosity and velocity especially at low velocities^ 

In order to fit successfully the total VDF from the 
present data we consider the following function that includes 



3 Here it is interesting to note that the early-type and the late- 
type populations are well fitted by Schechter luminosity functions 
both separately and in combination although they are well sepa- 
rated photometrically. 



<j)tot{<r)da = 0* 



+ £ 



x exp 



r(»//3) 

T(a'//3) 
13 da 



(£)" 



r(a//3) a 



where <f>* is the total number density and e represents the 
fraction of the number density due to the second (correc- 
tion) term. Fig. [2] shows that the total VDFs are well fitted 
by equation ((8j). For the fitting process there arises a weak 
degeneracy between parameters a and e. To break the de- 
generacy we fix the value of a to that for the late-type VDF 
since late-type galaxies dominate in the low velocity disper- 
sion limit. The fitted parameters of the total VDFs can be 
found in Table [T] 



3 STRONG LENSING STATISTICS: DATA 
AND MODEL 

3.1 Data 

For the purpose of this work we need galactic-scale strong 
lens samples that satisfy well-defined observational selection 
criteria. Such statistical lens samples that have been pub- 
lished so far include tho se from the radio-s e lected CLASS 
and PANELS surveys (|Mvers et all 120031 ; iBrowne et ail 
120031 ; IWinn et al l l200ll) a nd the optically-selected HST 
Snaps hot dMaoz et al.lll993l ) and SDSS Quasar Lens Search 
(QLS: llnada et al.ll2008l ) surveys. However, the SDSS QLS 
lens systems are limited to image separations greater than 
1 arc second while the other surveys include subarcsecond 
systems down to ~ 0".3. Moreover, the SDSS QLS lens 
systems are from relatively low-redshift (0.6 < z < 2.2) 
sources while the other surveys include high-redshift sources. 
For these reasons we use only the HST Snapshot and the 
CLASS /PANELS lens samples to study the evolutions of the 
VDFs. The HST Snapshot data include 4 lens systems out of 
506 quasars with z > 1. The HST Snapshot lens systems are 
summarised in Table [2] and the dat a for the entire sources 
can be found in lMao7et all l| 19931 ). We use the selection 
functi on for the Snapshot survey described bv lMaoz fc Rbd 
(1 19931) and th e sourc e number- magnitude relation described 
bv lKochanekl lll996l) to calculate magnification biases and 
cross sections. The CLASS data include 22 lens systems out 
of 16521 radio sources. However, only a subsample of 8958 
sources including 13 lens systems, referred to as the CLASS 
statistical sample, satisfies well -defined selection criteria (see 
IChaell2003l ;lB rownc et al.ll2003l ) so as to allow reliable calcu- 
lation of cross sections and magnification biases. We use the 
C LASS s electio n function and source properties described 
by IChaei (120031 ) except that we use the number-flu x density 
relation slightly updated bv lMcKean et all (|2007al ). The rest 
of 9 lens systems from the CLASS and additional 4 lens sys- 
tems from the PANELS (the southern sky counterpart of 
the CLASS) can be used for relative lensing probabilities 
(see IChadl2005l ). The CLASS and PANELS lens systems 
can be found in Table [2] 

Notice that in Table [2] both the image separation A8 
and twice the average radial separation (2Ri) are given for 
each lens system. Ad is the maximum possible separation 
amongst the image components. This is the value adopted 
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Table 2 . Summary of Stron gly-lensed Sources. The properties of the stron gly lensed s ystem s fro m the Snapshot op tical survey and the 
CLASS l|Browne et al]|2003h and PANELS radio surveys are revised from IKochanekl jl996l l and IChaej [|2003l . I2005T ). The first 13 radio 
sources (B0218+ 357 through B2319+051) are members of the well-defined CLASS statistical sample of 8958 radio sources ferowne et al.l 
2003; Chac 2003). The CLASS statistical sample along with the Snapshot sample are used fully for all their observed properties while 
the rest of the sources are used only for their relative image separation probabilities. Notations are the following: mg - B magnitude; /s - 
flux density at 5 GHz; z a - source redshift; Z; - lens redshift; AO - the maximum separation between any pair of the images; (2Ri) - twice 
the average separation of the images from the primary lensing galaxy (or the image centre if the lens is not observed) ; N- lm - number of 
images; E/L - ear ly/late- type galaxy; CL - lensing cluster/group; G2 - secondary (satelli te) lensing galax y. Re feren ces are the following: 
1 - th e CA STLES website (h ttp:/ /cfa-www.harvard .edu/ castles/); 2 -IChael ||2003|. [20051) ; 3 -IKochanekl Jl996t) ; 4 -ISiemiginowska, et al.l 
lll998ll: 5 -lYork et al l ll2005h: 6 - Biggs et al l d2003h: 7 -iRusin et al.l l|200ll);*8~ISurDi fc Blandfordl d2003h: 9 -ICohn et al.l 1 200 J): 10 - 
iMcKean et al.l ll2007bl); 11 -IChae et al l feOOll) ; 12 - iPhillips et al.l J2000T I; 13 - Iwinn et all 1 1|2003|) ; 14 - iLehar et al.l J2000h : 15 - iKing et al.l 



l|l997l l: 16 - iMore et all l|2008h 
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" Assumed early-type fo r their large inferred velocity dispersions of cr Ri 270 km s — 1 for H1413+117 and <r « 230 km s~ x for B1152+199. 
HFassnacht et all lll999l) suggests z a = 1.28 based only on a single broad emission line. This value for z a implies an uncomfortably large 
SIS velocity dispersion of a m 384 km s — 1 for Gl. However, lHamana et aL (2005) reports a measured central stellar velocity dispersion 
of c — 213 ± 23 km s — . An alternative value for z s would be 4.3 assuming the emission line is a Lyman alpha line. We take the latter 

value for this work. 



m previous studies (e.g. lChaell200il2005l L (2ft) is two times 
the average of the separations of the image components from 
the primary lensing galaxy position or the image centre if 
the lensing galaxy position is not known. If the lens were an 
SIS, the two values would be the same. However, for many 
systems the two values are somewhat different, (2Ri) usually 
being somewhat larger than A9. The systems with relatively 
large differences tend to be quadruply imaged systems with 
highly asymmetrical morphologies. The average values for all 
30 lens systems are A# = 1".38 and (2Ri) = 1".45 respec- 
tively while those for 9 lens systems that are only quadruply 



imaged are A8 — 1".45 and (2Ri) = 1".63 respectively. If 
we exclude 6 multiple-lens systems, the average separations 
for the rest 24 systems are AO = l".17 and (2Ri) = l".21 
respectively. Although the differences between AO and (2Ri) 
are not large, we consider the two possibilities in order to 
be as precise as possible. The question is then which of the 
two choices would correspond more accurately to the ve- 
locity dispersion of the primary lensing galaxy. Because the 
velocity dispersion is directly related to the Einstein radius 
Re ( i.e. the critical radius) via a lens model (see, e.g.. IChael 
2003), the answer may be found by comparing A8 and {2Ri) 
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with 2Re- The Einstein radius Re for a lens can only be re- 
liably determined from a detailed modelling of the system. 
This is particula rly true for highly asymm etrical systems 
with large shears. ICongdon fc Keetonl (J2005) model some of 
highly asymmetrical quadruple lens systems through a lens 
model with the isothermal radial profile and general angu- 
lar structur es. The highly asymmetri cal systems that are 
analysed bv lCongdon fc Keetonl (120051 ) include B0712+472, 
B1422+231 and B2045+265. As shown in Tabled for these 
systems the values for AO and (2 Rj) are quite differen t . If we 
take the Einstein radii found bv lCongdon fc Keetonl (120051 ) 
for these systems, then 2Re are between AO and (2Ri) but 
much closer to (2Ri). This means that (2Ri) is probably the 
better choice for this work. Hence the results from this work 
are based on (2Ri) rather than AO. 



Another important issue regarding the proper interpre- 
tation of the lens data is how to treat the multiple-lens sys- 
tems. The operational definition for a multiple-lens system 
to be adopted in this work is a system in which the cen- 
tres of distinct multiple lensing galaxies are found within 
the Einstein radius of the lens. For such a multiple-lens sys- 
tem care must be taken in interpreting the image separa- 
tion. If secondary galaxies are much less massive than the 
primary galaxy and/or lie well outside the Einstein ring, 
then secondary galaxies may contribute little to the Ein- 
stein mass so that the image separation may well corre- 
spond to the velocity dispersion of the primary galaxy up 
to small errors. Clearly, this appears to be the case at least 
for B1152+199 and B1030+074. On the other hand, if the 
lensing galaxies are of comparable mass and/or well within 
the Einstein radius as in interacting galaxies or a chance 
alignment (in projection) of galaxies at the same redshift 
(as in a group) or at different redshifts, then the image sep- 
aration cannot match well the velocity dispersion of the pri- 
mary galaxy. This appears to be the case for the majority 
of the multiple-lens systems. Table [3] summarises the Ein- 
stein radii of the lensing galaxies in the multiple-lens sys- 
tems. Indeed, the difference between the Einstein radius of 
the total lens potential approximated by (Ri) and that of 
the primary lensing galaxy R^' is quite large. This means 
that taking the face values of (2Ri) of the multiple-lens sys- 
tems would make a significant error in a statistical analysis 
of strong lensing. Simply ignoring the image separations of 
the multiple-lens systems would not be a perfect solution 
either because the primary galaxies of the multiple-lens sys- 
tems might be biased compared with 'single' lensing galax- 
ies. The average value of ' for the 6 multiple-lens systems 
in Table [3] is 0".79 while that for the rest of 24 systems is 
« (2Ri)/2 = 0".61. If B2108+213 were excluded assum- 
ing that the large separation is assisted by a group/cluster, 



then the average value of R^' would be 0".60. However, the 
measured velocity dispersion of Gl in B2108+213 appears to 
be consistent with the large separation w ithout any signifi - 
cant contribution from the group/cluster (|More et al.ll2008T ) 
implying that B2108+213 may well be a valid data point. 
Therefore, for the multiple-lens systems we use 2R^ instead 
of (2Ri) in our analyses. In doing so, we treat multiple-lens 
systems as if they were single-lens systems but with cor- 
rected image separations. 



3.2 Model 

Suppose that there is a population of cosmological sources 
described by a luminosity function or a number-flux den- 
sity relation N Zs (> /„) (number of sources with flux density 
greater than /„ in a certain comoving volume). Then, for 
a source with observed redshift z s and flux density /„ the 
probability that it is multiply-imaged (with a certain image 
multiplicity) by intervening galaxies is, crudely speaking, 

probability = distance to the source 

x number density of galaxies 

x multiple-imaging cross section 

x magnification bias (9) 

and can be formulated as follows in the order of the four 
factors: 





cdt 
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d(A0) 



F 

Jo 
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^ ct (A0,aO 



dN z ,(> f v /fi) dN Zs (> /„) 



dfu 



(10) 



where t is cosmic time (i.e. proper time of a comoving 
galaxy), (j>(a,z) is the velocity dispersion function, s(a, z) 
is the cross section of a galaxy with velocity dispersion a 
at redshift z, and A# m i n and ^i m i n are respectively the mini- 
mum allowed values of the separation and the magnification. 
In equation (|10[) . V a {AQ,n) is the probability density dis- 
tribution of image separation (AO) and magnification (/z) 
within the multiple-imaging region for a galaxy of velocity 
dispersion a. Va{A0, n) depends on the lens model. For the 
simple case of the singular isothermal sphere (SIS) it is given 
by 



V„{A6,ti) = 5{A0 
with 

A6>sis = 8tt 



A0si; 



, 2/i m i n 



Ta(z, Z s ) 



(!)' 



(11) 



(12) 



7\4(0, Z s , 

where c is the speed of light in vacuum and ta{z\,Z"i) is 
the angular-diameter distance between redshifts Zi and 22 
in units of the Hubble radius. 

We use the si ngular isoth ermal ellipsoid (SIE) lens 
model described bv lChael (|2003h to calculate lensing proba- 
bilities for a population of galaxies given by a VDF (equa- 
tion^. We assume a mean projected mass density axis ratio 
of 0.7 ignoring the scatter. The expressions for absolute lens- 
ing probabilities , i.e. t hose i ncluding magnification biases, 
can be found in IChael l|2003l ) while relati v e pro babilities of 
image separations can be found in Chac (2005). The total 
likelihood for lensing due to a certain population of galaxies 
is given by 



ln£ t o 



(13) 



ln£ 

opt + ln£ rad , 

where £ op t is the likelihood due to optically selected sources 
given by 

r (opt) 



U 



ln£ opt = ]T ln[l-p (opt) (fc)]+ Y, ln<5p ( ° pt) (Z), (14) 
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Table 3. Multiple-lens systems 



System 


Gl's ER 
R W 


2< 


G2's ER 


G3's ER 


(Ri) - 


Comments 


Reference 


B1152+199 


w 0.8 








fa 


G2 ignored 




B1359+154 


0.36 


0.72 


0.29 


0.29 


0.50 




Rusin et al. f 20011 


B1608+656 


0.53 


1.06 


0.29 




0.64 




Koopmans ct al. (2003) 


B2045+265 


1.06 


2.12 


0.09 




0.11 




McKean et al. (2007b) 


B2114+022 


0.89 


1.78 


0.46 




0.42 


two-plane lensing 


Chae et al. (2001) 


J0134-0931 


0.17 


0.34 


0.17 




0.24 




Keeton fc Winn (2003) 


B1030+074 


as 0.8 








fa 


G2 ignored 




B2108+213 


1.74 


3.48 


0.51 




0.55 


cluster /group ignored 


More et al. (2008) 



Notes: R 



(n) 



— n-th lensing galaxy's (Gn's) Einstein radius (ER); (Ri) - the average separation of the images from the primary galaxy 

(Gl). All the values are given in arcseconds. 



and £ ra d is the likelihood due to radio-selected sources given 
by 



COMPARING THE LOCAL VDFS WITH 
INTERMEDIATE REDSHIFT (0.3 < Z < 1) 
VDFS 



ln£. 



in[i-p (rad) (fc)]+ J2 inV rad) (o 
\ 



+ 



E ln ^ls ad) 0') 

3=1 



(15) 



Then a "x 2 " is defined by 



-2 In C to 



(16) 



Here p(k) is the integrated multiple-imaging probability, 
8p(l) is the differential lensing probability of the specific sep- 
aration, lens and source redshifts, and 5pis(j) is the relative 
image separation probability. 

For equation (|14[) the Snapshot data are used. The 
lensed systems are the 4 systems of the first part in ta- 
ble M and A^ opt) = 502. For equation (HU) the CLASS and 
PANELS data are used. For the first parenthesis in equa- 
tion (|15p the lensed systems are the 13 systems of the 2nd 
part of Table[2]and ATy d ^ = 8945. For the second parenthe- 
sis in equation (|15p the lensed systems are the 13 systems 
of the 3rd part of Tabled 

We calculate lensing statistics for three kinds of galaxy 
population, i.e. the early-type, the late-type and the total 
population. For the population of a certain morphological- 
type we must include only the lensed systems of the 
same morphological-type lensing galaxies. However, for sev- 
eral lensed systems the lensing galaxy types have not yet 
been determined by observations nor could be judged oth- 
erwise. These systems include Q1208+1011, B0128+437, 
B0739+366 and B1555+375. These systems all have rela- 
tively small image separations. Each of these systems may 
well be either early-type or late-type. Hence we include these 
systems for the statistics of both the early-type and the late- 
type populations but with their contributions to the \ 2 mul- 
tiplied by a penalisation factor of 0.5. 



For the source redshift distributions of the CLASS and the 
Snapshot the likely redshift of the lensing galaxy lies in 
0-3 < z < 1 (see Table [2}. Because of this, if we constrain 
a non-evolving VDF using the lens data the constrained 
VDF would correspond to an averaged VDF over the red- 
shift range or a VDF at z ~ 0.65 assuming the evolution 
is smooth in z. Comparing the constrained VDF of a cer- 
tain galaxy population with the corresponding local VDF 
may then reveal the essential features of the evolution of 
the VDF. 



We constrain the VDF of the form given by equation ((2)1 
for both the type-specific VDFs and the total VDF. The lo- 
cal total VDF is well-fitted by the six-parameter function 
of equation ((8]) rather than the function of equation ((2j) (see 
§2.4). However, the number statistics of the strong lens sam- 
ple is not strong enough to warrant the six-parameter func- 
tion. When all four parameters of the VDF (equation [2]| 
are allowed to vary, the confidence regions in the parameter 
space are broad because of parameter degeneracies. How- 
ever, the VDF itself is constrained well within plausible 
ranges of the parameters. Hence we may fix the low-velocity 
end slope a to the local value to break the parameter de- 
generacy without significantly altering the possible range of 
the VDF. 



Fig. [3] shows the 68% and 95% confidence limits (CLs) 
of the constrained VDF for the total and the type-specific 
galaxy populations. In Fig. [3] the local VDFs are overplot- 
ted and compared with the constrained VDFs. We notice 
the following trends for the evolutions of the VDFs. First, 
the number density of highest velocity dispersion galaxies 
appears to evolve most as can be seen from both the total 
VDF and the early-type VDF. Second, the lower velocity 
dispersion part of the early-type VDF does not show any 
signature of evolution. Finally, the late-type VDF evolves 
little overall. 




Figure 3. The strong lensing constrained intermediate-redshift (0.3 < z < 1) VDFs are compared with the local VDFs (Fig. [2] 
and table [TJ. These simple comparisons reveal the trends of the evolutions without detailed modelling. The green and red regions are 
respectively the 68% and the 95% confidence limits. 



5 CONSTRAINTS ON THE EVOLUTIONS OF 
THE VELOCITY DISPERSION FUNCTIONS: 
PARAMETRIC MODELLING 

We study the evolution of the VDF through a parametric 
approach. We parametrise the evolution of the number den- 
sity and the velocity dispersion as 

0» = <^1O F *; a*(z) = a *10 (C3/4)z , (17) 
and the evolution of the shape of the VDF as 

a(z) = ao (i + ^YT^) ; £(*)=#> { i + k i 3 YT~z) ' (18) 

Here z is the cosmological redshift and z/(l + z) — 1 — (o/ao) 
(where a and ao are the cosmological scale factors at red- 
shift z and zero respectively). Parameters P and Q in equa- 
tion (I17p are parametrised such that they have the same 
sensitivity to the strong lensing optical depth (the abso- 
lute lensing probability). The positivity of a and /3 in equa- 
tion (|18|) requires that k a > —1 and kp > — 1. Notice that 
parameters Q, k a and kp have sensitivities to the image 
separation distributions as well as the lensing optical depth 
while parameter P has sensitivity only to the lensing optical 
depth. 

The present lens sample is limited to the image sepa- 
ration Ad > 0".3. Consequently, the data cannot constrain 
the shape of the VDF for a < 95 km s _1 Q Hence the con- 

4 For a lens redshift of 0.6 and a source redshift of 2. A6 = 0".3 
corresponds to a as 95 km s -1 for the assumed cosmology. 



straints from this study are valid only for a > 95 km s _1 
and outside this range the results must be regarded as ex- 
trapolations. Because of this limitation of the lens data and 
the small number statistics parameter k a is ill-constrained 
by the data. Moreover, for a given shape of the VDF there 
are parameter degeneracies as pointed out in §4. Hence it is 
necessary to impose a constraint on the parameters. We may 
consider k a — or k a = hp as a default choice. The former 
choice corresponds to a constant shape for the low velocity 
end slope while the latter corresponds to a constant faint- 
end slope since ql = a(z)/f3(z) — 1 (equation[5]). In practise 
these two choices make little difference for the present data. 
We adopt the latter constraint because the intermediate- 
redshift VDFs constrained in §4 appear to be steeper both 
at the low velocity and the high velocity ends. 

5.1 Early- type galaxies 

Fig.Hshows the 68% and 95% CLs in the P-Q plane for two 
local early-type VDFs. The fitted parameter va lues can be 
found in Tableg] The early-type VDFs (both the lChoi etHI 
120071 ) measured E0 VDF and the Monte Carlo generated El 
VDF) evolve little overall. The overall little evolution of the 
early-type VDF is consistent with previous constraints from 
stron g lensing statistics (e.g. IChae fc Mad 120031 ; lOfek et"al] 
l2003t) . 

However, more detailed nature of the evolution as a 
function of velocity dispersion can be found in Fig. [S] in 
which the local VDFs are compared with the corresponding 
projected VDFs at z = 1. The main features of Fig. [5] may 
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Table 4. Evolution parameters for the velocity dispersion func- 
tions 
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Figure 4. The 68% and the 95% confidence limits in the evolu- 
tionary parameters of the early-type velocity dispersion function. 
The integrated number density evolves as 10 Fz while the velocity 
dispersion parameter a* oc lO^ 2 / 4 . The parameter for the evolu- 
tion of the shape is not shown here but can be found in table [4] 



be summarised as follows. First, the number densities of 
early-type galaxies of large velocity dispersions (a > 230 — 
250 km s~ ) have not only changed significantly from 2 = 1 
to z — but by an increasingly larger factor for a larger 
velocity dispersion. The dependence of the number density 
evolution on velocity dispersion appears to be real because a 
constant shape is clearly disfavoured by data. In particular, 
the number density of the largest velocity dispersion galaxies 
(er > 300 km s _1 ) has increased by a large factor (> 3) since 
z — 1. Second, for early- type galaxies of typical velocity 
dispersions and lower (a < 230 — 250 km s _1 ) there is no 
statistically significant change in number density from z = 1 
to z = 0. This explains why lens statistics gives the "right" 
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Figure 5. The early-type VDFs at z = 1 are compared with the 
local counterparts. The z = 1 VDFs are the functions projected 
using the constrained evolution parameters (table [4jl . The green 
and red regions show the 68% and 95% confidence limits due 
to the uncertainties in the constrained parameters. These figures 
clearly show the differential number density evolution for early- 
type galaxies of high velocity dispersion in which the evolution is 
greater at a higher velocity dispersion. 



cosmology assuming no evolution in the number density and 
the shape (see §6.2). 

Finally, we must bear in mind that the above results 
for the early-type population suffer from systematic errors 
arising from the lack of the knowledge of the lensing galaxy 
morphological type for several lens systems (see §3). The 
results for the late-type population shown below are prone 
to the same problem. However, the results for the evolution 
of the total VDF (§4.3) are free from the problem. 



5.2 Late- type galaxies 

The evolutionary behaviours of the late-type VDF can be 
found in Figs.[6]&[7] Overall, the late-type VDF is consistent 
with no evolution at the 68% confidence level. Furthermore, 
in contrast to the early-type VDF there is no statistically 
significant differential evolution for the late-type VDF. In 
particular, the number density of the most massive late- 
type galaxies does not change significantly between z — 1 
and z = 0. 
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Figure 6. The same as Fig. [4] but for the late-type population 
of galaxies. 
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Figure 7. The same as Fig. [5] but for the late- type population 
of galaxies. 
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Figure 8. The same as Fig. [4] but for the total population of 
galaxies. 



5.3 All galaxies 

We study the evolution of all galaxies as a whole using the 
total VDFs AO and Al in Table [1] One advantage of con- 
straining the evolution of the total VDF is that it is not af- 
fected by the uncertainties in the identification of the lensing 
galaxy type for several systems. However, parametric mod- 
elling of the total VDF is less straightforward because the 
total VDF has more parameters (equation [5} than the type- 
specific VDF (equation [2} . Namely, there arise more severe 
degeneracies amongst the evolutions of the parameters. To 
break the parameter degeneracies we choose to fix the pa- 
rameters of the correction term in equation i.e. a' and 
e. Fig. [8] shows the 68% and 95% CLs in the P-Q plane. 
The fitted parameter values can be found in Table [3] Notice 
that the fitted values of the evolution parameters are quite 
different from those of the early-type VDF because of the 
differences in the VDF functions as well as the data used. 

In Fig. [9] the total VDFs at z = are compared with 
the projected total VDFs at z — 1. The inferred evolution- 
ary trends of the total VDFs are consistent with those of the 
type-specific VDFs. There are only relatively minor quan- 
titative differences between the total VDF and the type- 
specific VDFs in the detailed nature of the evolution. Most 
noticeably, the number density evolution starts to become 
important at a « 200 km s _1 from the inferred evolution of 
the total VDF (whereas it starts at a ~ 250 km s _1 for the 
early-type VDF) and the degree of evolution for large ve- 
locity dispersion galaxies is greater for the total VDF than 
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Figure 9. The same as Fig. [5] but for the total population of 
galaxies. 



the early-type VDF. These minor quantitative differences 
are due to the differences both in the functional form of 
the VDF and the data used between the total VDF and 
the early-type VDF. The differences may reflect a possible 
range of errors. Despite the minor quantitative differences 
in the evolution between the total VDF and the early-type 
VDF, the overall trends of the evolution for large velocity 
dispersion galaxies are essentially the same. 



6 DISCUSSION 

6.1 Possible sources of systematic errors 

6.1.1 (The interpretations of) the strong lensing data 

The results presented in §5 are based on the most likely 
interpretation of the strong lensing data as judged by the 
author. However, the present data do not exclude alterna- 
tive interpretations of the strong lensing systems. We assess 
here the greatest possible effects on galaxy evolution of al- 
ternative interpretations of the strong lensing data. 

First, it is possible that the source redshift for 
B2045+265 is z a = 1.28 and accordingly the velocity dis- 
persion of the lensing galaxy is a « 384 km s _1 . We judge 
that this interpretation is unlikely based on the following 
arguments. First of all, the measu red velocity dispersion of 
the lensing galaxy as reported by lHamana et all (|2005h is 
a = 213 ± 23 km s _1 which is much lower than the im- 
plied SIS velocity dispersion. If the measurement is correct 



and the source redshift is z s = 1.28, this system would be 
an extremely rare case in which the SIS velocity dispersion 
is much larger than the measured stellar velocity disper- 
sion (ustars/crsis = 0.55 ± 0.06). For the Sloan Lens ACS 
Survey (SLACS) sample of 53 earl y- type lens galaxies with 
170 km s _1 < friars < 400 km s _1 [Bolton et all i|2008l ) find 
0.7 < (Jstars/o-sis < 1.3 with a mean value close to unity. Sec- 
ondly, according to lHamana et al" I (|2005l ) the lensing galaxy 
in B2045+265 lies in the fundamental plane with the mea- 
sured velocity dispersion implying that it is a normal ellipti- 
cal galaxy with the 'correct' velocity dispersion. Thirdly, the 
source redshift z s — 1.28 for B2045+265 would give a much 
worse x 2 for the model fit (A^ 2 > 12 for the evolution model 
of the total VDF). This means that the source redshift z s = 
1.28 for B2045+265 is not in harmony with all the other lens 
systems. Finally, the suggestio n of z s = 1.28 is base d only 
on a single broad emission line dFassnacht et al . 1999). This 
is far from a secure redshift from the observational point 
of view alone. Nevertheless, if we co nsider z s = 1. 2 8 and 
so take a — 384 km s" 1 ignoring the lHamana et all (|2005h 
measurement of a = 213 ± 23 km s" 1 for B2045+265, the 
evolutions of the early-type VDF and the total VDF would 
change significantly at the high end of velocity dispersion. 
In Fig. [10] the total VDF AO at z — is compared with its 
projected VDF at z = 1 with z s = 1.28 for B2045+265. The 
ratio 4>z=i/4>z=o decreases up to a ~ 300 km s~ but then 
turns upward so that it approaches near unity within the 
68% CL (although it is not shown) as a reaches 400 km s _1 . 
Namely, in this model the most massive early-type galaxies 
are already in place at z = 1 and there is little change in 
their number density since then. Overall the number den- 
sity of early-type galaxies as a function of velocity disper- 
sion is consistent with no evolution at 95% confidence be- 
tween z = 1 and z = 0. This is in stark contrast to the 
results presented in §5 without z s = 1.28 for B2045+265 
that require the greater evolution at the larger velocity dis- 
persion for a > 200 km s _1 . Given this potentially large sys- 
tematic error, it is important to measure more reliably the 
source redshift and c onfirm the stellar velocity dispersion of 
lHamana et ail (|2005l ) for B2045+265. 

Second, the environments of the lensing galaxies may 
cause biases in inferring the velocity dispersions of the lens- 
ing galaxies from the observed image separations. We cor- 
rected the effects of the secondary (and tertiary) lensing 
galaxies within the image regions based on the results of 
detailed lens modelling found in the literature (Table |3J). 
However, we ignored the possible effects of groups or clus- 
ters surrounding (and galaxies nearby) the image regions. It 
is thus possible that the velocity dispersions of the lensing 
galaxies were overestimated for some fraction of the lenses. 
There have been no systematic observations on the envi- 
ronments of the CLASS or the Snapshot lenses. However, 
a recent study of the environments of the SLACS lenses 
finds that typical contributions from the environments of 
the 'normal' lenses (i.e. those without very close compan- 
ions) are less than a few per cent in the lensing mass density 
l|Treu et alj|2008l) . For other several lens systems (including 
some CLASS and Snapshot lenses) with identified surround- 
ing groups/clusters recent studies find that the environments 
contribute typical ly less than 10% to the convergence of 
the lens potential dMomcheva et"a l. 2006; Auger et al. I l2007l ; 
iFassnacht et alj|2008l )7 Since the lens surface density is pro- 
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portional to the velocity dispersion squared, neglection of 
the environmental contributions would cause errors in the 
velocity dispersion at most about 3%. Theoretical simula- 
tions bv lKeeton fc Zabludofj (120041 ) predict somewhat larger- 
error of 6% for the velocity dispersion. To assess a maximum 
possible systematic error arsing from neglecting the lens en- 
vironments we consider lowering the observed image sepa- 
ration by 25% for the nine lens systems with a group (or 
groups) along the sight line (see table [2j| assuming that the 
implied velocity dispersion is biased by 5%. The evolution of 
the total VDF AO for this case can be found in Fig. [TO] The 
best-fit for this case gives a better x 2 (^X 2 = —2.5) com- 
pared with the case that the surrounding groups contribute 
nothing to the image separations. However, the evolutionary 
behaviour of the VDF is little changed. 

6.1.2 Statistical properties of sources 

For the statistical analysis of strong lensing the properties of 
the source population are of crucial importance. The statis- 
tical properties of the source population influence the lensing 
probabilities through its redshift distribution and number- 
magnitude (or number-flux density) relation or through an 
evolving luminosity function collectively. Consequently, any 
error in the source properties can bias the results based on 
strong lensing statistics. For the present work, we use two 
distinct populations of sources, namely optical quasars and 
flat-spectrum radio sources (CLASS sources). For the optical 
quasar population the source redshifts are all known and the 
number-magnitude relation is also known empirically. For 
the flat-spectrum radio source population the source prop- 
erties are less certain. For most of the radio sources the red- 
shifts are not known and the inferred redshift distribution 
for the source population was largely based on 27 measured 
redshifts for a subsample of sources with 5 GHz flux den - 
sities between 25 mjy and 50 mjy (|Marlow et al.| [2000). 
Since lensed sources can come from flux densities below 
25 mjy, it has been assumed for the previous analyses of the 
CLASS lensing that the redshift distribution is unchanged 
for sources below 25 mjy. Hence it is possible that the 
present work suffers from a systematic error in the CLASS 
parent source redshift distribution. However, IMcKean et al.l 
(2007a) reports (a preliminary result) that the mean red- 
shift for a sample of flat-spectrum sources with 5 GHz flux 
densities between 5 mjy and 15 mjy i s consistent w i th tha t 
for sources of higher flux densities by iMarlow et al. | l|2000h . 
The number-flux density relation of the flat-spectrum radio 
sources is well measured above 30 mjy at 5 GHz because it 
is based on a large number (~ 10 4 ) of source counts. It is less 
certain below 30 mjy because s uch a large number of s ources 
are not available. Nevertheless, IMcKean et al.l (l2007aT ) gives 
a measurement of the number-flux density relation below 
30 mjy based on a sample of 117 sources dow n to 5 mjy 
and c onfirms a previous estimate reported by IChae et al.l 
l|2002h . 

As discussed above there is no observational indica- 
tion that the source properties adopted in this work are 
grossly biased. However, the uncertainties in the CLASS 
source properties clearly require further observational stud- 
ies. It is thus of interest to see whether the two independent 
statistical lens samples, namely the Snapshot sample and 
the CLASS statistical sample, give agreeing results. The 
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Figure 10. Possible systematic errors in the evolution of the 
VDF. The green solid line is the best-fit for the AO VDF as 
shown in Fig. [9] The red dashed line is the fit assuming that 
the source redshift z s = 1.28 is adopted and interpreted to give 
a Ri 384 km s _1 for the lensing galaxy in B2045+265. The red 
dash-dot line is the fit assuming that the surrounding groups 
contributed significantly to the image separations and the cross 
sections for several systems (see the texts in §6.1.1). The red dot- 
ted and dash-dot-dot-dot lines are the fits respectively including 
only the CLASS sample and only the Snapshot sample for the 
absolute lensing probability. 



two independent results can provide a cross check of each 
other. Unfortunately, the Snapshot sample contains only 
four lenses so that it has a weak statistical power in con- 
straining the evolution of the VDF. In particular, the Snap- 
shot data alone cannot constrain the shape of the VDF at 
all. We thus include the image separation distribution (ISD) 
of all 26 radio lenses for the Snapshot fit. Notice that the 
ISD of the radio lenses does not depend on the statistical 
properties of the CLASS parent source population so that 
the inclusion of the radio ISD does not bias the result of the 
Snapshot fit. In Fig. [10] the radio and the Snapshot best-fit 
results are compared. Although not shown in the figure the 
68% CLs of the two results overlap so that they are consis- 
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tent with each other. However, the best-fit results indicate 
that the Snapshot sample requires more evolution overall 
but less evolution at the high mass end than the CLASS 
sample. 



6.1.3 Galaxy mass profiles 

The results of this work are based on the assumption that 
the galactic mass distribution follows the isothermal distri- 
bution in the inner region of the galaxy probed by strong 
lensing (within about the effective radius). This assump- 
tion naturally leads to the equality of the observed stellar 
velocity dispersion and the SIE model velocity dispersion. 
Conversely, if the SIE assumption is systematically violated 
the lens statistics cannot reliably constrain the observed lo- 
cal stellar velocity dispersion. Of course, not all galaxies will 
closely follow the isothermal mass profile. However, the de- 
viations of individual galaxies from the isothermality would 
not matter much statistically as long as the average profile 
is close to the isothermal. 

Systematic errors from the isothermality assumption 
can arise in the following two contexts. First, if the mass 
profile varies systematically depending on the velocity dis- 
persion, e.g. a higher velocity dispersion galaxy follows a 
shallower profile on average (or vice versa), then the results 
of this work can be systematically biased. The solution to 
this potential problem appears to be provided in part by the 
SLACS study of early-type galaxies through a combination 
of lens modelling and st ellar velocity dispe rsion measure- 
ment of lensing galaxies. iBolton et al.l (|2008h finds through 
an analysis of 53 massive early-type galaxies that the mass 
profiles of early-type galaxies scatter around the isother- 
mal profile but the average profile is close to the isother- 
mal profile independent of the v elocity dispers i on be tween 
175 km s" 1 and 400 km s" 1 . The IBolton et all (|200cf ) find- 
ing applies to the upper velocity dispersion part to the peak 
velocity dispersion (~ 160 km s _1 ) of the early-type VDF. 
For the early-type galaxies of the peak velocity dispersion 
and lower, it remains unanswered whether the isothermal- 
ity assumption is systematically unbiased. However, down 
to the velocity dispersion limit (a ~ 95 km s _1 ) sensitive 
to the strong lensing data (remember that this study is lim- 
ited by the image separation range) most early-type galaxies 
are likely to be normal elliptical/lenticular galaxies with lit- 
tle contamination of dwarf ellipticals/lenticulars (see, e.g., 
Ide Riicke et~aill2005l ). Hence it is likely that the isothermal- 
ity assumption remain not so bad for low velocity dispersion 
galaxies within the range probed by this work. For the late- 
type galaxies it is expected that the mass profile becomes 
systematically shallower in the inner region (< i?Jf|) as the 
luminosity decreases for sub-L« galaxies accordin g to the ro- 
tation curves of spiral galaxies (|Persic et al.lll996h . However, 
for the range of galaxies probed by the strong lensing data 
(circular rotation speed V c > 134 km s~ ) the deviations are 
not likely to be too large. The effects of the neglected errors 
can be understood as follows. The assumed constant rotation 
curve for a galaxy whose rotation speed actually declines to- 
ward the galactic centre will inflate the lensing mass and so 
the lensing cross section. This inflated lensing cross section 
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will then lead an underestimate of the number density so as 
to keep the observed lensing occurrences. Consequently, the 
evolution of the number density for the late-type popula- 
tion can be overestimated. We postpone the quantification 
of this systematic error to a future work although it may 
not be warranted because of the small number statistics of 
late-type lenses in the present data. 

Second, if the average mass profile of galaxies deviates 
systematically from the isothermal as z increases, then the 
results on the galaxy number density evolution from this 
work can be biased. The total (luminous plus dark) mass 
profiles of galaxies at intermediate and higher redshifts can 
be best constrained by modelling strong lens systems pro- 
vided with reliable and stringent observational constraints. 
The best studied examples in strong lensing appear to sug- 
gest that the avera ge total mass prof i le of g alaxies is isother- 
mal up to z ~ ll I lKoopmans et al. I (|2006r ) carry out a sys- 
tematic modelling of SLACS lenses and find that the mass 
profiles of gravitational lens galaxies are close to isothermal 
rega rdless of the lens redshi fts between z = 1 and z = 0. 
The iKoopmans et all l|2006T ) sample contains 6 lenses for 
0-5 < z < 1, so that the isothermality may be valid at least 
up to z ~ 0.7. Individual extensive modelling of lenses at 
z ~ 1 also find isothermal(- like) total mass distributions 
of galaxies. ISaha et al.l l|2007l ) find a nearly isothermal pro- 
file from a non-parametric m odelling of MG0414+054 at 
z = 0.96. IWinn et all (|2004T l model J1632-0033 at z ~ 1 
and find a tight constraint of j3 = 1.91 ± 0.02 for a mass 
density of p(r) oc r - ' 3 (/3 = 2 correspon ding to the isother- 
mal) using an observed central image. iTreu fc Koopmana 
(2002) combines galactic dynamics and lens modelling to 
find P = 2.0 ± O.l(stat) ± O.l(syst) for MG 2016+112 at 
z — 1.004. Moreover, there exist many modelling exam- 
ples in which the authors assume the isothermal profile for 
the lensing galaxies at z ~ 1 and reproduce the observed 
image properties although some systems appear to require 
substructures. To sum up, strong lensing studies support 
the isothermal profile up to z ~ 1 without any indication 
of evolution of the total mass profile in redshift. The to- 
tal mass profiles can also be constrained by other dynam- 
ical pr obes although few studies hav e been carried out at 
z > 1. Ivan der Wei fc van der Marell (|2008l ) present veloc- 
ity dispersion profiles of early-type galaxies up to ~ 4 kpc 
(0".5) at z ~ 1 which are consistent with isothermal profiles. 

Studies based on surface photometry find varying re- 
sults on the evolution of the galaxy mass profiles. Many 
recent studies suggest that high redshift (z > 1) galax- 
ies may be more compact in their light and stellar mass 
density distributions than the local cou n terparts (e.g. 
Truiillo et ai1l2007l; Ivan der Wei et~ai1 120081; Icimatti et all 
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particular, Bez anson et al 



2008; Bezanson et al 



2009). 



In 

(|2009l ) find that the stellar 
density profiles at z ~ 2.3 are more steep ly declining be- 
yond r = 1 kpc than the local counterparts. Bez anson et al.l 
(2009) also note that the stellar density within 1 kpc is 2-3 
times higher at z ~ 2.3 than at z ~ 0. This structural evo- 



6 It is now well established that lensing galaxies at intermedi- 
ate redshifts z ~ 0.3 — 0.7 have isothermal (or close to isother- 
mal) mass profiles. The question here is whether the isothermality 
would hold up to z = 1. 
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lution may b e explained by galaxy evolution through mi- 
nor mergers ilNaab et all 120091 120071 ). On the other hand, 
iHopkins et all (|2009l ) find similar densities for the observed 
physical radiu s rang e both at z > 2 and at z ~ 0. 
Hopk ins et all (|2009l ) further note that there may or may 
not be differences in the densities outside the central regions 
between z > 2 and z ~ 0. The ambiguity outside the cen- 
tral regions arises primarily from the difficult y in observing 
the low surface brightness parts. Interestingly. ISaracco et aD 
(2009) find that early-type galaxies at 1 < z < 2 may be 
divided into two groups, one group showing no size evolu- 
tion and the other group implying a significant size evolu- 
tion compared with local galaxies. These studies may imply 
that at least some galaxies undergo size evolutions possi- 
bly including stellar mass profile evolutions. Complicating 
the issue, dark mass profiles may also evolve as a result of 
(responding to) stellar mass growth and dark mass growth 
via merging. Unfortunately, the baryon-involved evolution 
of dark mass profile is a poorly understood part of galaxy 
evolution. Consequently, the surface photometry data alone 
cannot settle the issue of the total mass profile evolution 
even if the stellar mass profile evolution is real. 

As described above direct probes of the galaxy mass 
profiles give support for the isothermality up to z = 1. How- 
ever, the evidence is not broad and strong enough to rule 
out the variation from the isothermality in z and/or a (or 
mass) . Hence galaxy mass profiles remain a source of uncer- 
tainty. If the true mass profiles deviate systematically as a 
function of z and/or a, then the number density evolution 
in z as a function of a will be biased. Namely, the variation 
in mass profile leads to the variation in the cross section and 
the magnification bias resulting in a different number den- 
sity evolution. Roughly speaking, if the true mass profile is 
steeper/more concentrated (shallower/less concentrated) at 
2 = f than the isothermal profile, then the number density 
at z = 1 should be lower (higher) than that derived here 
based on the isothermal profile. It is unclear at present how 
important the effects turn out to be quantitatively. A real- 
istic ray-tracing simulation would eventually be required to 
quantify the effects robustly. 



6.1.4 Local velocity dispersion functions 

The local velocity dispersion function provides a benchmark 
for the study of galaxy evolution just as the local luminos- 
ity function is. Hence we cannot overemphasise the impor- 
tance of reliable local VDFs for this work. For this reason 
we have considered n ot only the measured early-type VDF 
bv IChoi et al l (E007h but also the Monte Carlo simulated 
VDFs for both the early-type and the late-type populations 
based on the galaxy counts from the SDSS and the intrinsic 
correlations between luminosity and internal velocities. 

Despite the range of the VDFs considered in this work 
one may still consider the following systematic errors for 
the local VDFs. First of all, the galaxy counts for the type- 
specific galaxy populations remain somewhat un certain due 
to the error in the galaxy classifications. The [Ch oi et al.l 
(2007) classification matches relati vely well the morp hologi- 
cal classification. However, even the lChoi et al.l d2007l) classi- 
ficatio n has an error of about 10% according to lPark fc Choil 
(|2005h . 

Second, the early-type VDF is less reliable at lower ve- 



locity dispersions (a < 150 km s" 1 ) because of the incom- 
pleteness in the volume and magn itude limited samples of 
galaxies used bv lChoi et alJ |2C)07i) to obtain the early-type 
VDF and the intrinsic correlation between luminosity and 
velocity dispersion. Indeed, a Monte Carlo simulated SDSS 
galaxy sample makes a difference in the early-type VDF at 
velocity dispersions < 125 km s _1 (see Fig. [2}. However, the 
present strong lensing data are sensitive only to velocity dis- 
persions > 95 km s _1 so that the difference in the VDF for 
a < 125 km s _1 makes little difference in the results on the 
evolution of the VDF. 

Finally, the late-type VDF is relatively less reliable. We 
turned the circular rotation speed to the velocity disper- 
sion assuming the SIS profile for all late-type galaxies, but 
the observed mass profile deviates systematically from the 
isothermal toward the galactic centre. In this respect the 
late-type VDF may not be well-defined and it may be even- 
tually necessary to use directly the circular velocity function 
for the statistics of strong lensing through a realistic lens 
model of a late-type (spiral) galaxy. Such a realistic mod- 
elling will become very important in the future when the 
number statistics of late-type lenses becomes large. 



6.2 Comparison with previous results in lensing 
statistics 

The importance of galaxy evolution in strong lensing statis- 
tics was rec ognised and appl ied to a handful of l enses for the 
first time bv lMad <|l99ll) and lMao fe Kochanekl l| 19941 ). How- 
ever, strong lensing statistics has usually been used to con- 
strain the cosmological constant (and dark energy) and ap- 
pears to confirm the concordance cosmology according to re- 
cent results assuming the ea rly-type population evolves little 
between z = 1 and now (e.g.lChae et al.ll2002l ; lMitchell et ajj 
120051 ; IChaell2007l ; lOguri et al.ll2008l ). The assumption of no- 
evolution of the early-type population appeared also to be 
supported by strong lensing statistics assuming the con- 
cordance cosmological model acc ording to an analy sis of 
the CLASS statistical sample by IChae fc Maol (|2003l ) and 
a lens-redshift test bv lOfek et atl i|2003l ). Notice that these 
studies assumed a constant shape of t he the inferred veloc- 
ity dispersion function . In other words, IChae fc Mad (|2003T ) 
and lOfek et al.l l|2003l ) considered only the evolution of the 
number densit y and the characteris t ic vel ocity dispersion. 
More recently, ICapelo fc Nataraianl l|2007h carried out the 
lens-redshift test more extensively and studied in particular- 
possible systema tic errors based on a sampl e of 70 galaxy 
lenses. However, ICapelo fc Nataraianl l|2007l ) also assumed 
non-evolving shape for the VDF and obtained results con- 
sistent with earlier results. To sum up, all previous results on 
galaxy evolution from strong lensing statistics were based on 
the assumption of the non-evolving shape for the VDF and 
appeared to be consistent with no evolution of early-type 
galaxies. 

In this work we allow the shape of the VDF to evolve 
and find that strong lensing statistics actually requires a 
differential evolution of the VDF (see Figs. [S]and|5J). The 
differential evolution is intriguing in that the number den- 
sity of galaxies changes little between z — 1 and now be- 
low a ~ 200 km s~ but at higher velocity dispersions the 
number density evolution becomes significant with an in- 
creasingly larger factor for a higher velocity dispersion (see 
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Fig. [9)l . We discuss the implications of this differential evo- 
lution below in §6.4. Here we just point out that this be- 
haviour of differential evolution may explain why previous 
lens statistics resulted in little evolution assuming a constant 
shape. For a constant shape of the VDF its evolution may 
be realised only through the vertical and horizontal shifts, 
namely the variations of the integrated number density and 
the characteristic velocity dispersion. The shifts of the VDF 
are likely to be determined by the behaviour of the number 
density at the lensing peak velocity dispersion, i.e. the ve- 
locity dispersion for which the lensing probability (oc 0cr 4 ) 
is maximal. However, as shown in Fig.|5]the number density 
of galaxies evolves little at intermediate velocity dispersions. 
Consequently, the VDF appeared to evolve little by previous 
lensing statistics. 

6.3 Comparison with theoretical predictions and 
galaxy surveys: a qualitative analysis 

In the ACDM hierarchical structure formation picture the 
dark halo mass function (DMF) evolves in cosmic time as 
a consequence of hierarch ical merging (e.g. IWhite fc Reesl 
1 19781 ; lLacev fc Cold ll993T ). Accordingly, the stellar mass 
function (SMF) and the velocity dispersion function (VDF) 
of galaxies also evolve. However, baryon physics complicates 
the evolutions of the SMF and the VDF making it non-trivial 
to compare the evolutions of the DMF, the SMF and the 
VDF one another. Conversely, careful analyses of the coevo- 
lution of the DMF, the SMF and the VDF may reve al key 
insigh ts into galaxy formation and evolution processes ijChad 
2010). Here we compare the evolutionary trends of the VDF 
from this work against the DMF from N-body simulations 
and the SMFs predicted by recent semi-analytical models 
of galaxy formation and obtained by cosmological surveys 
of galaxies. This is meant to be an inter-comparison of the 
qualitative features of three functions as suggested by cur- 
rent theoretical and observational studies. A comprehensive 
and quantitative analysis of the coevolution of three func- 
tions is given in IChad (|2010h where evolutions of functions 
(j(Af v j r ) and a(M s t a rs) are derived and their implications 
are discussed in the context of current cosmological obser- 
vations. 

A generic prediction of the ACDM theory is the differ- 
ential evolution of the DMF in which the number density 
evolution of a larger M v ir (virial mass) halo is greater, e.g. 
since z — 1 (see Fig. [TTJ. As can be noticed easily, this evo- 
lutionary trend is strikingly similar to that of the total VDF 
(Fig. [9]). In fact, this qualitative similarity is expected once 
we assume that a is on average a monotonically increasing 
function of M v i r at each epoch (a plausible assumption) . The 
quantitative detail of the evolution of u(M v i r ) has implica- 
tions on the structural evolution of the galaxy-halo system 
|Chadl2010T l. 

Unlike the DMF the evolution of the SMF is not well 
predicted by theoretical models because of complex galaxy 
formation processes. Nevertheless, it appears that recent 
semi-analytic models (SAMs) predict a hierarchical evolu- 
tion of the SMF that shows evolutionary patterns similar 
to those of the DMF and the SMF (see Fig. [12]). How- 
ever, the 'observed' SMFs from galaxy surveys show different 
evolutionary patterns (Fig. I12[l . In particular, many galaxy 
surveys suggest stellar mass-downsizing (apparently anti- 
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Figure 11. The evolution of the dark halo mass function (DMF) 
between z = and z = 1 under the curren t ACDM par adigm. The 
results are based on the code provided bv lReed et ail (l2007t) . No- 
tice that the DMF evolution is differential for M vir > 10 13 h _1 M Q 
as is the VDF evolution for a > 200km s _1 (see Fig.[9]l- 
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20081 : iMarchesini et al 
2009h . although there are results that 
larly support stellar mass-downsizing (e.g 
iFaber et al.ll2007l : lBrown et al.l[2008l ; lllbert et al.ll2009l ). The 
disagreement between SAM predictions and galaxy survey 
results is mor e apparent for lower ste llar mass galaxies. Some 
authors (e.g. iFontanot et al.l 120091 ) suggest that this dis- 
agreement reflects some additional physical processes miss- 
ing in the current SAMs. However, assuming that a is on 
average a monotonically increasing function of M stars at 
each epoch, we may expect a hierarchical evolution of the 
SMF. In this view, the VDF evolution is qualitatively more 
in line with the current SAM predictions than the down- 
sizing SMFs from galaxy surveys. However, if downsizing 
SMFs from galaxy surveys should turn out to be true, a 
deeper rethinking of cr(M ata rs) may be required. A detailed 
quantitative analysis of the evolution o f cr(M v j r ) ba sed on 
currently available data is carried out in IChad (|2010l ). 



7 SUMMARY AND CONCLUSIONS 

In this work we have constrained the evolutionary be- 
haviours of the velocity dispersion functions of galaxies 
based on the local measured and Monte Carlo realised 
VDFs and the statistics of strongly-lensed radio-selected and 
optically-selected sources. This work is also based on the as- 
sumption that the total (luminous plus dark) mass profile is 
isothermal in the optical region for ^ z ^ 1 in line with 
strong lens modelling results. The constrained evolutions of 
the VDFs can be characterised as follows. 

• The (comoving) number density of massive galaxies 
with a > 200 km s~ , which are mostly early-type galaxies, 
evolves differentially from z = 1 to z — in the way that 
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the evolution is greater at a higher velocity dispersion. The 
most massive galaxies of a ~ 400 km s _1 are much rarer 
at z = 1 than the present epoch with its number density 
probably less than 30% of the local density. 

• The number density of intermediate-mass galaxies with 
95 km s" 1 < a < 200 km s" 1 are nearly constant between 
z = 1 and z — regardless of the morphological type. 

• The velocity dispersion function of late-type galaxies 
transformed from the Monte Carlo realised circular velocity 
function is consistent with no evolution between z = 1 and 
z = either in the overall number density or its shape. 

Comparing the constrained VDFs with the ACDM 
DMF and the SMFs observed from galaxy surveys as well as 
predicted from semi-analytical models of galaxy formation, 
we find the following. 

• The VDF evolution is in line with the evolution of the 
ACDM DMF, i.e., the hierarchical build-up of mass struc- 
tures over cosmic time. 

• The evolutionary pattern of the VDF is similar to the 
SMF evolution predicted by recent semi-analytical models 
of galaxy formation but qualitatively different from the stel- 
lar mass-downsizing evolution suggested by galaxy surveys. 
Further investigations are required to clarify this apparent 
conflict. 

In conclusion we have investigated in detail the evo- 
lutionary behaviours of the VDFs for the first time based 
on local SDSS data and statistical strong lensing data and 
found that the evolution pattern is strikingly similar to that 
of the ACDM DMF. Also, there is a promising agreement be- 
tween the VDF evolution and the SMF evolution predicted 
by SAMs although observed stellar mass-downsizing evolu- 
tion needs yet to be clarified. Larger well-defined statistical 
samples of strong lensing from future observation tools like 
the SKA will be invaluable laboratories for decoding galaxy 
formation and evolution processes. 
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Figure 12. (a) Stellar mass functions (SMFs) of galaxies at z = 
0. Blue solid, dashed, dash-dott ed and dott e d cur v es are respec- 
tively t he observed SMFs fr o m lllbert et al.l 1I2OO9I) . iPanter et al.l 
|(2003), IConroy fe Wechslerl ll2009l l and iBell et al.l ll 2003l l. No- 
tice that the Iconrov fc Wechslerl d2009l > SMF is a composite 
of sev eral galaxy surv e y res ults. Black solid curve is a median 
of the IFontanot et alj feoogf l three semi-analytic model (SAM) 
predictions wh il e bla ck dashed curve is the prediction of the 
ICattaneo et al.1 1120081 ) SAM. (b) SMFs at z = 1. The colours 
and styles of the li nes refer to the s ame r eferences as in (a). Black 
dotted curve is the lStringer et aL 1 1|2009|) SAM prediction. Notice 
the discrepancy between the observed SMFs and the predicted 
SMFs for the lower part of the stellar mass range, (c) SAM pre- 
dicted SMFs at 2 = and 2 = 1 are compared. Black curves 
refer to 2 = while red curves refer to 2 = 1. Solid and dashed 
curv es are respective l y the predictions by IFontanot et al.l (2009) 
and ICattaneo et al.1 J2OO8). Notice that these predicted SMFs 
show hierarchical differential evolutions in line with the VDF and 
the DMF (see Figs. [9] and [TTJ. (d) Observed SMFs at 2 = and 
2 = 1 are compared. Blue curves refer to 2 = while red curves 
refer to 2 = 1. Solid and dash- dotte d curves are respectively th e 
SMFs from lllbert et al.l j2009h and IConrqy fc Wechslerl l|2009h . 
The composite SMF bv lConrov fc Wechslerl (120091 7 shows an anti- 
hierarchical evolution while the COSMOS SMF by lllbert et all 
II2009T ) does not show such a differential evolution. 



